Ballistic transport in graphene beyond linear response 
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The process of coherent creation of particle - hole excitations by an electric field in graphene 
is quantitatively described beyond linear response. We calculate the evolution of current density, 
number of pairs and energy in ballistic regime for electric field E using the tight binding model. 
While for ballistic flight times smaller than tni oc E~^^^ current is linear in E and independent of 
time, for larger ballistic times the current increases after tni as J cc E^^^t and finally at yet larger 
times {t > tB cc E~^) Bloch oscillations set in. It is shown that the number of pairs follows the 
2D generalization of the Schwinger's creation rate n cc E^l'^ only on certain time segments with a 
prefactor different from that obtained using the asymptotic formula. 

PACS numbers: 81.05.Uw 73.20.Mf 73. 23. Ad 



INTRODUCTION 

It became increasingly evident that electronic mobil- 
ity in graphene is extremely large exceeding that in best 
semiconductor 2D samples [l[. Since the system is so 
clean the transport becomes ballistic, especially in sus- 
pended graphene samples 0] , so that interactions of elec- 
trons with phonons, ripplons, disorder and among them- 
selves can be neglected. Therefore there is a chance to 
observe various theoretically predicted exotic phenomena 
like nonlinear response and Bloch oscillations 0]. 

Generally ballistic transport occurs due to two distinct 
phenomena. If a mobile charge carrier is available (as in 
an electron gas in a metal), electric field accelerates it, 
so that current increases linearly in time. In addition, an 
electric field can create new charge carriers (the process 
typically suppressed by energy gaps). A peculiarity of 
the ballistic transport in graphene with the Fermi level 
pinned right on the two Dirac points [J] (which happens 
naturally in suspended graphene samples [2|) is that there 
are no charge carriers present at all. The Fermi surface 
therefore shrinks to just two points. The carriers are 
created solely by an applied electric field like in the Zener 
tunneling effect in semiconductors but energy gap 
vanishes due to "ultra - relativistic" dispersion relation, 
e = |k| , where Vg ^ 10^™/ sec is the graphene velocity. 

The electron - hole pairs are created fast enough to 
make the current linear in electric field and constant 
in time, so that it looks like a Drude type linear re- 
sponse due to disorder rather than the ballistic accelera- 
tion J oc t of an electron gas with finite carrier density. 
The illusion of the "Ohmic" behaviour however cannot 
continue indefinitely in the absence of scatterers, and 
should eventually cross over to some sort of "accelera- 
tion" or even Bloch oscillations at large times. The be- 
haviour is expected to become nonlinear as function of 
electric field as indicated by the nonlinearity of the pair 



creation rate. It was shown long time ago in the 
context of particle physics, that the pair creation rate at 
asymptotically large times is proportional to E^^^ . 

Ambiguities in the application of the standard Kubo 
approach for the ultra - relativistic spectrum led us 
propose a dynamic approach to the tight binding model 
of graphene Within leading order in E (linear re- 
sponse) wc found that the DC conductivity is cr2 = f x 
rather the often cited value obtained both from 

Kubo formula Q and within the Landauer formalism [§] . 
In this note we solve the tight binding model for arbitrary 
constant electric field. The evolution of current den- 
sity demonstrates that the crossover from the "Ohmic" 
regime to the nonlinear one occurs at the experimentally 
achievable time scale tni oc E~^^^. Bloch oscillations are 
shown to set in on scale Ib oc E~^ much longer than tni 
for experimentally accessible electric fields. We discuss 
the relevance of the 2 D ge neralization of the Schwinger's 
creation rate formula [lO| to physics of graphene. 



TIGHT BINDING MODEL AND ITS EXACT 
SOLUTION 

Electrons in graphene are described sufficiently well for 
our purposes by the 2D tight binding model of nearest 
neighbor interactions on a honeycomb lattice [J. The 
Hamiltonian in momentum space is 
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with 7 ~ 2.7eV being the hopping energy and sum is 
over the Brillouin zone. Nearest neighbors are separated 
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by distance a = 3 A, 5 = 2 cos {ak^/ 2) and pseudospin 
index denotes two triangular sublattices. We consider 
the system in a constant and homogeneous electric field 
E along the y direction switched on at < = 0. It is de- 
scribed by the minimal substitution p = fik + - A with 
vector potential A = (0, —cEt) for t > 0. Since the cru- 
cial physical effect of the field is a coherent creation of 
electron - hole pairs, mostly near the two Dirac points, 
a convenient formalism to describe the pair creation is 
the "first quantized" formulation described in detail in 
0, [m ■ To consider the ballistic transport at zero tem- 
perature, r = dynamically, one starts at time t = 
from the zero field state in which all the negative en- 
ergy one - particle states, — |/ik| = — Ekj are occupied. 
The second quantized state evolving from it is uniquely 
characterized by the first quantized amplitude. 
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FIG. 1; The time evolution of the current density (3 bottom- 
curves) and the scaled pair creation rate (3 top curves) at 
relatively short times for various electric fields. 
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which is a "spinor" in the sublattice space. It obeys the 
matrix Schroedinger equation 
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It is a peculiar property of the tight binding matrix 
Eq.([T]) that solution for arbitrary ky can be reduced to 



that for ky ~ and has the Fourier series: 



V'k(i)= Yl Pmexp(-^<^); 
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where t ^ t — t^aky/S and ujf^ = + 3flm for frequency 
n = £/ {2V3ty) ; £ = E/Eq. The relevant microscopic 
time scale is tj = h/j and field £'0 = 7/ (ea). Recursion 
relations for the Fourier amplitudes Pm , 
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has two solutions p''*, s = ±1 with two Floquet frequencies 
[1^ Lo^ . The recursion is easily solved numerically and 
has the following convergent expansion in b in the whole 
relevant range, Q < b < 2, 
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with Wq = sfl+yt-y ^ -I- O^. It turns out that the two Flo- 
quet frequencies obey the relation obeying w"*" = 20— w^. 



again peculiar to graphene, as can be checked by both 
the perturbation theory, Eq.([7]) and numerical results. 
For experimentally accessible cases O << and the 



frequencies are just 
initial conditions 



Coefficients arc fixed by 



V'k (t = 0) = Uk 
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This solution is used to calculate evolution of current 
density, energy and the number of electron - hole pairs. 



TIME SCALE FOR OBSERVATION OF THE 
BLOCK OSCILLATIONS IN GRAPHENE 

Evolution of the current density during the ballistic 
" flight time" thai is the integral over Brillouin zone (mul- 
tiplied by factor 2 due to spin) 0: 



Jyit)^-2eJ2i'l (t) 
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The current density divided by electric field, a [t] = 
Jy (t) /E, is shown in Fig.l and 2 for various values of 
the dimensionless electric field £ in the range 2~* — 2~^. 

Fig. 1 in which evolution is shown up to ballistic time 
of 120t-y, demonstrates that after an initial fast increase 
on the microscopic time scale t-y (shown in more detail, 
using linear response, in Q), cr (i) approaches the uni- 
versal value (72 and settles there. Beyond linear response 
one does not expect the current density to hold up to 
this value indefinitely. In a ballistic system the energy 
initially increases, as follows from the Joule law. The 
total energy of electrons can be written in the first quan- 
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FIG. 2; The time evolution of the current density is shown 
for various fields up to 1000 



tized formalism as 

Utot (i) = 2 ^ it) Hpiji, (t) = 2 iij (t) \H\ ^ (t)) . 

(10) 

It can be shown using Eq.(|4]) that the power 



P{t) 
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= EJy (t) , 



is indeed proportional to current density. Since in the 
tight binding model electron's energy cannot exceed the 
upper band edge energy 87, hence at some time scale 
Ib the energy increase is reversed. The physics which 
takes over is that of the Bloch oscillations and is similar 
to that in ordinary materials, namely, electrons' energies 
are elevated by the electric field [3| due to the quasi- 
momentum shift. This feature is not related to the unique 
" rclativistic" feature of the graphenc spectrum. 

The current, shown in Fig. 2 for ballistic times up to 
lOOOt.^, indeed exhibits Bloch oscillations. It turns out 
that the current vanishes at points given exactly at mul- 
tiples oUb/^ with 
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being the period of the Bloch oscillations. The Bloch 
time is approximately the time required for the electric 
field to shift the momentum across the Brillouin zone 
Apy = eEtB ~ h/a. These times are very long for exper- 
imentally achieved fields, much longer than the ballistic 
flight time. One observes in Fig. 2 another peculiar fea- 
ture that (apart from the "rclativistic" initial constant 
segment) time dependence of a it) is similar for different 



electric fields. Indeed, if one plots JjyE versus tE^ all 
the curves nearly coincide. Moreover 



J{t)^ri<J2E]I^E^I^sx^[ — 



tB 



(13) 



is an excellent fit. 

For a sample of submicron dimensions, L ~ 0.5/im, 
W = 1.5/im, the ballistic time can be estimated as 



tbal = L/Vg 



2.3 • with 



For cur- 



rent as large as /max = "m-A the electric field is i?max = 
= lO'^V/m corresponding to £ = 10^"^ (voltage 
in such case would be quite large Vmax = 5V). The 
first maximum of the Bloch oscillation will be seen at 
flight time of ^3/4 = 3.6 • lO^t-y, which is of the same 
order as thai- If one uses a value of the current typi- 
cal to transport measurements / — bOfiA, the electric 
field is just i? = 5 • lO^V/m corresponding £ = 5 • 10"^, 
<b/4 = 7.2 • 10H~f >> tbal and is therefore out of reach. 



THE CROSSOVER FROM LINEAR TO 
NONLINEAR REGIMES 

In Fig.l one clearly observes a remarkable feature: 
there is a much smaller crossover time tni, after which 
the conductivity rises linearly with time above the con- 
stant "universal" value a2- 



J (t) = (72 



The crossover time is 
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It becomes the same as the ballistic time t^ai = 2.3-10"^t^, 
mentioned above, for relatively weak fields E = W^V/m 
corresponding to £ = 10~^. Therefore some of the trans- 
port measurements performed might be influenced by the 
physics beyond linear response. 

A qualitative picture of this resistivity without dissipa- 
tion is as follows. The electric field creates electron - hole 
excitations mostly in the vicinity of the Dirac points in 
which electrons behave as massless relativistic fermions 
with the graphene velocity Vg playing a role of velocity of 
light. For such particles the absolute value of the velocity 
is Vg and cannot be altered by the electric field and is not 
related to the wave vector k. On the other hand, the ori- 
entation of the velocity is influenced by the applied field. 
The electric current is ev, thus depending on orientation, 
so that its projection on the field direction y is increased 
by the field. The energy of the system (calculated in 
a way similar to the current) is increasing continuously 
if no channel for dissipation is included. Therefore the 
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" Ohmic" conductivity originates in creation of pairs near 
the Dirac points with an additional contribution due to 
the alignment of the particles' motion with the field's 
direction. At times of order Ib his process exhausts it- 
self due to the following processes. Electrons gain mo- 
mentum from the electric field and leave eventually the 
neighborhoods of the Dirac points. They are no longer 
ultra-relativistic and are described by (positive or nega- 
tive) effective mass [1] and the more customary physics 
takes over. 

The crossover to the nonlinear regime can be detected 
from within the perturbation theory in electric field. In- 
deed wc found that the E'^ correction to conductivity is 



J (t) /E = a2 
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The correction therefore becomes as large as the leading 
order for t — 2.1tj/£^^^ ~ tni- To gain more insight 
into the nature of the crossover to nonlinear response we 
calculated also evolution of the energy and number of 
electron - hole pairs during the ballistic flight. 



SCHWINGER'S PAIR CREATION FORMULA 
AND GRAPHENE 

The states in the conduction band for each momentum 
k in the Brillouin zone are described by a pscudospinor 
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orthogonal to defined in Eq.®. The amplitude of 
lifting of an electron into the conduction band is — 
(-0 (t) |wk) and consequently the density of pairs (factor 2 
for spin) reads, 



iVp(0 = 25^|Ak|^ = 2^ 



'^2 



(18) 



and the rate -^Np is shown in Fig.l as function of time. 
Its time dependence exhibits several time scales. At times 
smaller than expansion in electric field is applicable 
and the leading order result is: 
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This is analogous to "linear response" for current. Im- 
mediately after the switching on of electric field (times of 
order t^) the behaves as t^. For < t < tni the pair cre- 
ation rate per unit area rises linearly (with logarithmic 
corrections) , 



dt " 
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and is dominated by the neighborhood of the Dirac 
points. 

However it is clear from Fig.l that the expansion 
breaks down at tj-a, when the rate stabilizes approxi- 
mately at 



3/2 



and scales as the power E^^^. The rate continues to 
rise in a series of small jumps till Bloch oscillations set 
in. At that stage (actually at about t_B/4) number of 
electrons elevated into the conduction band becomes of 
order one, consistent with Eg. p^ . Then it oscillates. 
The power E^^'^ is, up to a constant, the same as the 
rate of the vacuum breakdown due to the pair produc- 
tion calculated beyond perturbation theory by Schwinger 
in the context of particle physics (when generalized to the 
2-t-l dimensions and zero fermion mass [1, [l3|)- This is 
not surprising since the power i?"^/^ is dictated by di- 
mensionality assuming ultra - relativistic approximation 
is valid. However the physical meaning is somewhat dif- 
ferent. Wc have used here a definition of the pairs num- 
ber with respect to Fermi level of the system before the 
electric field is switched on (equivalently when an elec- 
trons are injected into a graphene sheet from a lead). 
This is different not only from the Schwinger's path in- 
tegral definition in which the Fermi level is "updated" 
along the work of electric field and from the definition 
proposed recently [lo| in connection with graphene. The 
asymptotics at very large times is not relevant for exper- 
imentally achievable ballistic times, so that the predicted 
relatively short plateau segments are more important. 



SUMMARY 

Ballistic transport in single graphene sheet near Dirac 
point was investigated using the dynamic approach be- 
yond linear response theory. We found that, while the 
observation of the Bloch oscillations is difficult, there ex- 
ists a novel time scale tni, see Eq.((T5]), of transition to a 
nonlinear regime which is within reach of current experi- 
mental techniques. The physics of the ballistic transport 
in graphene can be described as a succession of four time 
segments with different character. 

(i) At microscopic ballistic times t ^ t-y the current 
reacts fast to electric field and depends on microscopic 
details. 

(ii) The current density at zero temperature stays con- 
stant a2E for ballistic times < t < tni and physics 
is partially universal in the following sense. There are 
generally two contributions to the current. While one 
contribution is dominated by Dirac points, the other is 
related to the band structure. However the second contri- 
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bution vanishes due to symmetry properties of the Bril- 
louin zone, see ref. 

(iii) For t„i < t < ts the current density during the 
flight would rise above this value. It is dominated solely 
by the close vicinity of each of the two Dirac points. 
Perhaps the increase of conductivity might be at least 
partly responsible for the "missing tt" problem 0, ITst . 
namely that experimentally measured minimal conduc- 
tivity is higher than 172 even in suspended samples 0]. 

(iv) Finally at t ^ ts Bloch oscillations set in. The 
physics is again dominated again by the band structure, 
is "non - relativistic" and is not directly related to the 
Dirac points. 

It should be noted that in addition to limitations of 
the tight binding model used which ignores impurities, 
interactions, deviation of the chemical potential from the 
Dirac point and temperature beyond linear response such 
"relativistic" effects like the pair annihilation neglected. 
For very large electric fields the effects of radiation of 
energy into space (radiative friction) might in principle 
be observable and should be investigated. On the other 
hand influence of temperature and nonzero chemical po- 
tential in nonlinear regime are expected to be similar to 
those in linear response studied in Q. 

We are grateful to S. Gurvitz, E. Andrei, A. Mor- 
purgo, J. Goltzer, W.B. Jian for discussions. Work of 
BR and YK was supported by NSC of R.O.C. grant 
#972112M009048 and MOE ATU program. BR acknowl- 
edges the hospitality of the Applied Physics Department 



of AUCS; M.L. acknowledges the hospitality and support 
at Physics Department of NCTU. 



* Electronic address: |vortexbar@yaIioo.com| 

[1] S.V. Morozov et al, Phys. Rev. Lett, 100, 016602 (2008). 

[2] X. Du, I. Skachko, A. Barker and E. Y. Andrei, Nature 
Nanotechnology 3, 491 (2008). 

[3] J. H. Davis, "The physics of low dimensional semicon- 
ductors", Cambridge University Press (1998). 

[4] A. H. Castro Neto et al. Rev. Mod. Phys. 81, 109 (2009). 

[5] J. Schwinger, Phys. Rev. 82, 664 (1951). 

[6] K. Ziegler, Phys. Rev. B75, 233407 (2007). 

[7] M. Lewkowicz and B. Rosenstein, Phys. Rev. Lett, 102, 
106802 (2009). 

[8] E. Fradkin, Phys. Rev. B 33, 3257 (1986); PA. Lee, 
Phys. Rev. Lett, 71, 1887 (1993); V. P. Gusynin and S. 
G. Sharapov, Phys. Rev. B 73, 245411 (2006). 
[9] M.I. Katsnelson, Eur. Phys. J. B 51, 157 (2006); J. 
Tworzydlo et al, Phys. Rev. Lett. 96, 246802 (2006). 

[10] A. Casher, H. Neuberger and S. Nussinov, Phys. Rev. D 
20, 179 (1979); D. AUor, T. D. Cohen and D.A. McGady, 
Phys. Rev. D 78, 096009, (2008). 

[11] E.S. Fradkin, D.M. Gitman and S.M. Shvartsman, Quan- 
tum Electrodynamics with Unstable Vacuum (Springer- 
Verlag, Berlin 1991). 

[12] C. Chicone, "Ordinary differential equations with appli- 
cations". Springer - Verlag, New York (1999). 

[13] K.S. Novoselov et al. Nature 438, 197 (2005); Y. Zhang 
et al, Nature 438, 201 (2005). 



